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1. INTRODUCTION 


The reporting period is characterized by a change and reorganisation of the 
data distribution as well as of our data processing system available at the 
various institutions engaged in image processing in Zurich. 

Data reception and distribution switched from GSFC to the Telespazio S.A. and 
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its ground station in Fucino near Rome, Italy (announced with letter by NASA 
of January 14 th and hy Telespazio of April 14th, 1976). 

Numerous contacts took place hut up to now no- definite procedure for the data 
delivery could be arranged and no data was received. In addition to these or- 
ganisational problems the creation of a new format for image data files (as 
described in Quarterly Report No 5 ) caused quite some technical problems. ' 

The adaption of all software to the new situation, which is still under pro- 
gress, is wearisome and time consuming, hut necessary for further progress. 

Nevertheless significant scientific progress could be achieved and a new prd-jeet 
started during this period in the field of 

- snow mapping; 

- spectral reflectance measurements of snow on the ground; 

- attempt for automated comparison of multitemporal data for land use purposes 
by geometric correction of the data (new project). 


2. ACCOMPLISHMENTS 


2.1 Data processing 


Adaptation of the existing software to the new data format, uniform for all 
research groups engaged in image processing in Zurich (Zurich data format), the 
testing of the various preprocessing and classification systems as well as the 
installation and setting in operation of the new PDP-ll/ 40-RAMTEK-display system 
available now for interactive image analysis kept us busy for most of our time. 
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2.2 Snow mapping 


2 

A larger *e3t area of about ,2*700 km in Grisons (Pratigau - Schanfigg - Land- 
schaft Davos - Albula - Lenzerheide - Rheintal) Was classified based on our 
.software package (as described in Quarterly Report l/2, 1975) for LANDSAT frame 
E-2090-09234 of April 22nd, 1975* 

For the supervised digital classification with the PENN STATE ORSER SYSTEM 
(PPD-CLASS) 19 different categories (learning groups) combined into 5 main 
categories (Fig. l) were set up. Particular interest was given to the melting 
zone (transition zone) with snow patches of various sizes intermixed with snow- 
free surfaces (vegetation, rocks, bare soil, roads etc.) to classify and map the 
position of the temporary snow line as exact as possible. The classification' 
could be achieved by using MSS channels 5 and 7 only (because no bright snovfree 
features were present in the test area) * 

To locate the temporary snow line as best as possible careful ground controls 
and aerial under flights during the LANDSAT overflight were undertaken. The data 
then was correlated with the satellite data. For the classification output the 
Optronics Photomation System was used, which allows a geometric correction and 
a selection of an optimal scale within certain limitations (see Quarterly Report 
No 3, 1976). With this method it is possible to locate the transition, zone bet- 
ween the compact snow cover and the totally snowfree areas with high accuracy 
for all different locations (relief, exposure, shadows etc.). The same test area 
was interpreted by analog methods on small-scale aerial photographs of the same 
date and the temporary snow line transferred onto a map 1:200' 000 using a ZOOM 
Transfer Scope. 

Fig. 2 shows the distribution of the 19 categories in the- feature plane of LANDSAT- 
channels 5 and 7 by means of ellipsis, where the center presents the mean and the 
two main half-axis the standard deviation of each category. 

Fig. 3 contains the classification- output of the total test area for the 3 main 
categories "snow - transition zone - snowfree areas" in a scale approx. 1:240' 000 
and Fig. 4 an enlarged section of the region of Davos in a scale approx. 1:64 '500. 
Both images were produced with the Optronics Photomation System. In Fig. 4 the 
transition zone in the valley bottoms as well a3 above the tree line is clearly 
indicated by its light gray tone. 



Fig* 5 is an: output with the line printer of a section of the same area as in 
Fig 4, using the same 5 main categories. The .rather "broad "band of the transition 
zone representing the "snow line" is deary recognisable. 

Quantitative .studies will be added on the areal extent of the boow cover and 
the transition" zone and on the altitudinal position of the temporary snow line 
for the total area as well as in its regional differentiation and its correlation 
with other features (exposure, slope angle, surface elements etc.). 


2.5 Spectral reflectance of snow measured with the EX0TECH-100 


Measurements of various snow surfaces were taken at the ground and from a heli- 
copter with the EXOTECH-lOO', At the same time the global radiance was measured^ 
or the reflectance was calibrated by using a reference panel. Prom these data 
the target reflectance was calculated to allow a direct comparison of the 
results. 

There is a clear correlation with the wavelength. For all snow types the target 
reflectance decreases with longer wavelengths between 0,5 - 1,1 jm and from 
MSS band 4 to 7 (Pig. 6). 

A general correlation between the target reflectance and the condition of the 
snow surface or with the snow parameters measured on the ground (depth, water 
equivalent, temperature, roughness of the surface etc.) could not be established. 

Therefore no definite conclusions can be made on the influence of these parameters 
on the target reflectance. On the other hand certain trends are obvious. A signi- 
ficant chance of the target reflectance occurs with the grain form(to be recog- 
nized best in channel 4), the snow depth (up to 50 cm), the wetness (to be 
recognized best with ratio 4/7 or 5/7), the pollution of the surface (channel 4) 
and partly with the roughness of the surface (channel 7). 

All these parameters should be taken into account for further measurements . For 
reflectance measurements from a helicopter additional variables have to be con- 
sidered such as atmospheric transmission, scattered light, relativ humidity etc. 
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A separation of the various snow types based on the target reflectance is 
possible to a certain extent. The fresh snow (sticky snow and slush snow) 
can be definitely separated* from the old snow (spring snow, brittle snow, 
brittle crust and solid crust) except for old dry snow (Fig. 6). A better 
• classification of snow types is only possible if additional snow parameters 
such as age, pollution, depth etc. are included, too. Fig. 6 shows the target 
reflectance for different snow types (mean value of measurements) in the- 4 
M3S-channels. 


2.4 Geometric correction of LANDSAT data for multitemporal analysis 

A new project was started to use multitemporal LANDSAT data in an automated 
classification procedure of land use features. The test area is the "Grosses 
Moos" near the Lake of Neuchdtel and the Lake of Bienne -in Western part of 
Switzerland. To exactly correlate and compare pixels of the same area from 
different orbits geometrical transformations are needed. To identify exact 
ground points small reflecting mirrors were constructed and set up during 
LANDSAT overflights, to get a clear mark on the image of these ground points. 
Difficulties occur because the exact orbit parameters of the satellite passes 
are not known previous to the overflight. We have been asking Telespazio S.A. 
for assistance in this respect but no information was received until today 
whilst the experience had to be carried out in June. Approx, orbit parameters 
were derived from past satellite passes by calculating a mean orbit. 



3. SIGNIFICANT RESULTS 


Using MSS channels 5 and. 7 and a supervised classification system with a PPD 
classification algorithm (Penn State Univ. ORSER System PPD-CLASS) it was 
possible to map the exact areal extent of the snow cover and of the transition 
zone with melting snow patches and snow free parts of various ' sizes over a 
large area under different aspects such as relief, exposure, shadows etc. 

A correlation of the data from ground control, aerial underflights and earth 
resources satellites provided a very accurate and detailed qualitative as well 
as quantitative interpretation of the melting procedure of snow in high moun- 
tain environment, allowing a comparison of the extent and position of the 
transition zone in its regional variation and its correlation to exposure, 
slope angle, surface features etc. 
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Fig. 1 ; Categories and main categories considered- for an exact location of the 
transition zone (after H.P. STIRNEKANN ) 
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Fig. 2 : Presentation of the selected 19 categories (numbeis as in Fig. l) 

in form of ellipsis for LANfSAT-channe Is 5 and 7 
(after HJP. STIRNEMAHN) 
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Fig. 4 : Digital mapping of snovr and transition zone for the region of Davos, 

Grisons; scale 1:64*500 (after H.P. STIRNEMAJIN) 







Fig* 5 : Output of PPD-classification of a section of Davos , Orisons (section 

included in Pig. 4 ) with MSS-channels 5 and 7 (after H.P, STIRNEMOf) 
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Pig. 6 : • Reflectance of various snow types in the 4 LARDSAT-KSS-channels 
measured with. EX0TECH-100 (after K. STABNZ) 
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1. INTRODUCTION 


The reporting period is characterized by a change and reorganisation of the 
data distribution as well as of our data processing system available at the 
various institutions engaged in image processing in Zurich. 

Data reception and distribution switched from GSFC to the Telespazio S.A. and 
its ground station in Fucino near Rome, Italy (announced with letter by NASA 
of January 14th and by Telespazio of April 14th, 1976). 

Numerous contacts took place but up to now no definite procedure for the data 
delivery could be arranged and no data was received. In addition to these or- 
ganisational problems the creation of a new format for image data files (as 
described. in Quarterly Report No j) caused quite some technical problems. ' 

The adaption of all software to the new situation, which is still under pro- 
gress, is wearisome and timeconsuming, but necessary for further progress. 

Nevertheless significant scientific progress could be achieved and a new project 
started during this period in the field of 

- snow mapping; 

- spectral reflectance measurements of snow on the ground; 

_ attempt for automated comparison of multitemporal data for land use purposes 
by geometric correction of the data (new project). 

2. ACCOMPLISHMENTS 


2.1 Data processing 


Adaptation of the existing software to the new data format, uniform for all 
research groups engaged in image processing in Zurich (Zurich data format) , the 
testing of" the various preprocessing and classification systems as well as the 
installation and setting in operation of the new PDP-ll/ 40-RAMTEK-display system 
available now for interactive image analysis kept us busy for most of our time. 
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2.2 Snow mapping 


£ 

A larger test area of about 2 '700 km in Grisons (Pratigau - Schanfigg - land- 
schaft Davos - Albula - lenzerheide - Rheintal) was classified based on our 
software package (as described in Quarterly Report l/2, 1975) for LANDS AT- frame 
£-2090-09234 of April 22nd, 1975. 

For the supervised digital classification with the PENN STATE ORSBR SYSTEM, 
(PPD-CLASS) 19 different categories (learning groups ) oombined into 5 main 
categories (Fig. l) were set up. Particular interest was given to the melting 
zone (transition zone) with snow patches of various sizes intermixed with snow- 
free surfaces (vegetation, rocks, bare soil, roads etc.) to classify and map the 
position of the temporary snow line as exact as possible. The classification 
could be achieved by using MSS channels 5 and 7 only(because no bright snowfree 
features were present in the test area) . 

To locate the .temporary snow line as best as possible careful ground controls 
and aerial underflights during the LANDSAT overflight were undertaken. The data 
then was correlated with the satellite data. For the classification output the 
Optronics Photomation System was used, which allows a geometric correction and 
a selection of an optimal scale within certain limitations (see Quarterly Report 
Mo 3, 1976). With this method it is possible to locate the transition zone bet- 
ween the .compact snow cover and the totally snowfree areas with high accuracy 
for all different locations (relief, exposure, shadows etc.). The same test area 
was interpreted by analog methods on small-scale aerial photographs of the same 
date and the temporary snow line transferred onto a map 1:200' 000 using a ZOOM 
Transfer Scope. 

Fig. 2 3 hows the distribution of the 19 categories in the feature plane of LAMDSAT- 
channels 5 and 7 by means of ellipsis, where the center presents the mean and the 
two main half-axis the standard deviation of each category. 

Fig. 3 contains the classification- output of the total test area for the 3 main 
categories "snow - transition zone - snowfree areas" in a scale approx. 1:240' 000 
and Fig. 4 an enlarged section of the region of Davos in a scale approx. 1:64’ 500. 
Both images were produced with the Optronics Photomation System. In Fig. 4 the 
transition zone in the valley bottoms as well as above the tree line is clearly 
indicated by, its light gray tone. 



Fig. 5 is an: -output with the line printer of a section of the same area as in 
Fig 4 , using the same 3 main categories. The rather broad band of the transition 
zone representing the "snow line" is deary recognizable. 

Quantitative studies will be added on the areal extent of the snow cover and 
the transition zone and on the altitudinal position of the temporary snow line 
for the total area as well as in its regional differentiation and its correlation 
with other features ( exposure 9 slope angle, surface elements etc.). 


2.3 Spectral reflectance of snow measured with the EXOTECH- 100 


Measurements of various snow surfaces were taken at the ground and from a heli- 
copter with the EX0TECH-100. At the same time the global radiance was measured^ 
or the reflectance was calibrated by using a reference panel. From these data 
the target reflectance was calculated to allow a direct comparison of the 
results. 

There is a clear correlation with the wavelength. For all snow types the target 
p§fi©et§fie§ d§er§3g§§ with longer wav§i§agths Q,S = 1,1 jm and £?e® 

MSS band 4 to 7 (Fig. 6). 

A general correlation between the target reflectance and the condition of the 
snow surface or with the snow parameters measured on the ground (depth, water 
equivalent, temperature, roughness of the surface etc.) could not be established. 

Therefore no definite conclusions can be made on the influence of these parameters 
on the target reflectance. On the other hand certain trends are obvious. A signi- 
ficant chance of the target reflectance occurs with the grain form(to be recog- 
nized best in channel 4)» the snow depth (up to 30 cm), the wetness (to be 
recognized best with ratio 4/7 or 5/7), the pollution of the surface (channel 4 ) 
and partly with. the roughness of the surface (channel 7). 

All these parameters should be taken into account for further measurements. For 
ref lectance , measurements from a helicopter additional variables have to be con- 
sidered such as atmospheric transmission, scattered light, relativ humidity etc. 



A separation of the various snow types based on the target reflectance is 
possible to a certain extent. The fresh snow (sticky snow and slush snow) 
can be definitely separated from the old snow (spring snow, brittle snow, 
brittle .crust and solid crust) except for old dry snow (Fig. 6). A better 
classification of snow types is only possible if additional snow parameters 
such as age, pollution, depth etc. are included, too. Fig. 6 shows the target 
reflectance for different snow types (mean value- of measurements) in the 4 
MSS-channels . 


2.4 Geometric correction of LANDSAT data for multitemporal analysis 


A new project was started to use multitemporal LANDSAT data in an automated 
classification procedure of land use features. The test area is the "Grosses 
Moos" near the Lake of HeuchAtel and the Lake of Bienne in Western part of 
Switzerland. To exactly correlate and compare pixels of the same area from 
different orbits geometrical transformations are needed. To identify exact 
ground points small reflecting- mirrors were constructed and set up during 
LANDSAT overflights* to get a elear mark on the image of these ground points. 
Difficulties occur because the exact orbit parameters of the satellite passes 
are not known previous to the overflight. We have been asking Telespazio S.A. 
for assistance in this respect but no information was received until today 
whilst the experience had 'to be carried out in June. Approx, orbit parameters 
were derived from past satellite passes by calculating a mean orbit. 
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3. SIGNIFICANT RESULTS 


Using MSS channels 5 and 7 and a supervised classification system with a PPL 
classification algorithm (Penn State Univ* ORSER System PPD-CLASS) it was 
possible to map the exact areal extent of the snow cover and of the transition 
zone with melting snow patches and snow free parts of various sizes over a 
large area under different aspects such as relief, exposure, shadows etc* 

A correlation of the data from ground control, aerial underflights and earth 
resources satellites provided a very accurate and detailed qualitative as well 
as quantitative interpretation of the melting procedure of snow in high moun- 
tain environment , allowing a comparison of the extent and position of the 
transition zone. in its regional variation and its correlation to exposure, 
slope angle, surface features etc. 
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Fig. 1 : Categories and main categories considered for an exact location of the 

transition zone (after H,P. STIRNEKANN) 
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Fig. 2 : Presentation of the selected. 19 categories (numbeis as in Fig. l) 

in form of ellipsis for IAN DS AT- ch anne la 5 and 7 
(after H.P. STIRNEMANN) 
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Fig. 4 : Digital mapping of snow and transition zone for the region of Davos, 

Grisons; scale 1:64’500 (after H.P. STIRHEMANN) 





Pig. 5 Output of PPD-classification of a section of Davos,' Orisons (section 
included in Pig. 4 ) with MSS-channels 5 and 7 (after H.P. STIRN&MM3T) 
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Fig. 6 : Reflectance of various snow types in the 4 LANDS A T-K5 3 - channel 3 

measured with EXOTECH-lOO (after K. STAENZ) 




